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We describe a calculation of the on-shell renormalization factors in QCD and QED at the three 
loop level. Explicit results for the fermion mass renormalization factor Z m and the on-shell fermion 
wave function renormalization constant Z2 are given. We find that at 0{a 3 a ) the wave function 
renormalization constant Z2 in QCD becomes gauge dependent also in the on-shell scheme, thereby 
disproving the "gauge-independence" conjecture based on an earlier two-loop result. As a byproduct, 
we derive an Oipti.) contribution to the anomalous dimension of the heavy quark field in HQET. 



I. INTRODUCTION 



Pcrturbative calculations in field theory are often instrumental in establishing accurate relations between theory and 
experiment. For example, in heavy quark physics, that now occupies the central place in the study of CP violation, 
perturbative effects play an important role alongside with non-pcrturbative effects. For heavy quarks, one separates 
soft and hard contributions through an expansion in the inverse quark mass; familiar examples being HQET (for a 
review see e.g. I]]) and NRQCD j2|. It often happens that the leading term in such an expansion is determined by a 
perturbative calculation in full QCD and for this reason perturbative calculations with heavy on-shell quarks in the 
initial and final state received significant attention in recent years. For many processes of interest we have witnessed 
a rapid computational progress and currently the "standard" level of accuracy is the next-to-next-to-leading or 0{ot 2 s ) 
order. Further progress in this direction will require three-loop computations. 

In QED, an interesting recent development was an application of dimcnsionally regularized NRQED to the calcula- 
tion of fundamental properties of simple atoms, like positronium and hydrogen. Here too, higher precision will require 
on-shell three loop calculations in many cases. 

The renormalizability of QCD and QED implies that any multi-loop calculation can only be made meaningful if the 
corresponding renormalization constants are knownQ; hence the three-loop computations require the knowledge of the 
three- loop renormalization constants in QCD. As is well known, the renormalization constants are renormalization 
scheme dependent. In some schemes, as e.g. in the MS scheme, they have been calculated up to the fourth order in 
perturbation theory. This is however insufficient if one is interested in the processes with heavy quarks being on-shell 
in the initial and final state - precisely as in QED, the renormalization constants in the on-shell scheme are required. 

If one adopts the on-shell renormalization scheme for the quarks in QCD, two specific renormalization constants 
need to be calculated. The mass renormalization constant Z m is defined by mo = Z m m, where mo is the bare and m 
is the pole mass. The wave function renormalization factor Z 2 is defined as ipo = \fZ2%\), where ■0o an d V 1 arc the bare 
and renormalized quark fields respectively. Both renormalization constants can be obtained from the expression for 
the heavy quark propagator close to the mass shell. The mass renormalization constant Z m is determined from the 
position of the pole of the heavy quark propagator, while Z 2 is determined from the residue. Because of the infrared 
catastrophe, the residue does not exist and the introduction of the infrared regulator becomes necessary. In what 
follows dimensional regularization is used to regulate both ultraviolet and infrared divergences which appear in the 
calculation. 
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: An immediate consequence of a more precise knowledge of just the on-shell renormalization factors is that it leads to an 
exact 3- loop result for the relation between the pole-quark-mass and the MS-quark-mass, which was recently presented in M . 
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Let us remark that the choice of the infrared regulator is not important for the pole mass, since it is known to 
be infrared finite and gauge invariant to all orders in perturbation theory. On the contrary, Z2 is not infrared finite 
already in the first non-trivial order in perturbative expansion. Hence, choosing the infrared regulator merely amounts 
to defining Z2 (in addition to usual scheme definition) and one may expect that the question of gauge invariance may 
then depend on that choice^. Indeed, that is exactly what happens. An interesting phenomenon occurs if dimensional 
regularization is applied to infrared divergences. In this situation, in QED, Z2 becomes gauge invariant to all orders 
in the coupling constant. This feature follows from dimcnsionally regulated Johnson- Zumino identity 



dlog£ 2 , a f d D k 



d£ °7 {2ir) D k A 

where £ is the gauge fixing parameter in a general covariant gauge. Note that the right hand side of the above equation 
is zero by virtue of the fact that scale-less integrals in dimensional regularization are defined to be zero. We will derive 
Eq.(^) in Section VI using the path integral formalism. 

Though no relation similar to Eq.(Q) has been derived in QCD, the use of dimensional regularization for infra-red 
divergences resulted in some interesting observations. The one-loop Z 2 in QCD is gauge invariant since it is a trivial 
generalization of the QED result. The two-loop contribution to Z2 is already sensitive to the non-abelian structure 
of QCD and there is no a priori reason to expect it to be gauge invariant. However, an explicit computation in Ref. 
produced a gauge invariant result. The authors of that reference then conjectured that this fact may be true to 
all orders in perturbation theory, though no supporting arguments have been given. 

The purpose of this paper is to present the three loop calculation of Z m and Z2 in dimensional regularization and 
hence to give a complete set of three loop renormalization constants that are specific to the on-shell scheme. The 
knowledge of Z m and Z2 with such an accuracy is an important step towards three loop calculations with heavy quarks 
close to the mass shell. The bulk of the paper is technical and attempts to provide some details about how the actual 
computation is done. Though the scale of this calculation makes it impossible to provide a detailed account of the 
algorithms used to obtain the result, we at least try to highlight some of its aspects which might be of importance for 
future related work. A peculiar feature of our result is that, in contrast to two first orders in perturbative expansion, 
the Q(a s ) contribution to Z2 turns out to be gauge dependent thus disproving the all-orders gauge independence 
conjecture of Ref. @. 

The paper is organized as follows. In the next Section we introduce the necessary notations and describe the 
calculation. In Section III we discuss in detail how new master integrals that appear in non-abelian theory are 
computed. In Section IV the results for QCD renormalization constants are summarized. In Section V the anomalous 
dimension of the heavy quark field in HQET is derived. In Section VI we give the three loop renormalization constants 
for QED. In the appendices a complete list of master integrals for the three loop on-shell two-point functions, as well 
as some other useful formulas, are presented. 

II. PRELIMINARIES 

We are going to calculate the heavy quark mass and the heavy quark wave function renormalization constants 
in the on-shell renormalization scheme. The bare quark mass mo and the bare quark field ipo are renormalized 
multiplicatively : 

too = Z m m, tpo = \fZ^,ip, (2) 

where m and ip stand for the pole mass and the renormalized quark field, respectively. Both renormalization constants 
can be derived by considering one-particle irreducible quark self-energy operator £(p, m). We discuss this derivation 
in some detail since we believe that discussion given in |Q is slightly over-complicated in that place. 

Because of the Lorentz invariance, the one-particle irreducible quark self energy operator to) can be parame- 
terized by two independent functions: 

E(p, to) = mSi(p 2 , to) + (p — to)£ 2 (J} 2 , m), (3) 



2 In what follows, we loosely use the terms "gauge invariant", "gauge independent" and "gauge parameter independent" as 
equivalent. 
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so that complete fermion propagator reads: 



p — m Q + h(p, mj 

Let us expand the self energy operator in formal Taylor series around p 2 = to 2 : 



(4) 



E(p,m)« to Ei(p 2 ,m)| 2 _ 2 + m-^ Ei(/,to)| 2 2 (p 2 - to 2 ) + (p - to) E 2 (p 2 ,to)| 



mT,i(p z ,m)\ p 2 =m 2 +(p-m) ( 2m 2 ^ Ei(p 2 , m)| p2=m2 + E 2 (p*,m)| ?r m _ 



<9p 2 



(•5) 



Substituting this expression into the fermion propagator Eq.(Q), we identify the position of the pole with the pole 
mass and the residue with the wave function rcnormalization constant. We obtain: 



Z m = 1 + Ei(p 2 ,TO)| p2=m2 , 
^ = 1 + 2to 2 Ae i(p 2 ,to) 



2 2 + T>2(p , m) 2 2 



(6) 



Both of these constants can be obtained easily once an expansion of E(p, to) close to the mass shell is available. 
Let us introduce the Minkowski vector Q with Q 2 = m 2 and consider a one-parameter family of external momenta 
p = Q(l + 1). The quark self-energy operator reads: 



E(p, m) = toEi(p 2 ,to) + (Q — to)E 2 (p 2 , to) + tQE 2 (p 2 , to). 



Taking the trace 



Tj = Tr 



Q + m 



Am 2 



and expanding in f up to C(i) we obtain: 



Ti = E!(p 2 



E(p,m) 



<9 



Ei(p 2 ,?n) +tT, 2 (p 2 ,m), 



2m dp^ ' m) lp 2 =^ 



+ Eofp , to) , 



C(t 2 ). 



(7) 



(8) 



(9) 



Comparing Eqs.(^) with Eq.(^), we see that expanding the self energy operator in Qt along with projecting on the 
relevant structure by taking the trace as shown in Eq.(|8|) delivers an expression one needs to reconstruct the on-shell 
mass and wave function renormalization constants. The pole mass in the above equation is calculated iteratively 
by using mass counterterms where appropriate, i.e. by calculating lower order diagrams with the appropriate mass 
counterterm insertion. 

Hence, both Z m and Zi can be extracted from the heavy quark self-energy operator evaluated on-shell. This 
task is rather straightforward in low orders of perturbation theory but it becomes increasingly difficult when one 
goes to higher orders. Both the number of diagrams one has to calculate and the complexity of integrals increases 
dramatically. In the present case we have about sixty diagrams to be calculated and the question about an efficient 
way to do the calculation becomes of tremendous importance. 

The most efficient way to evaluate those multiloop integrals is to utilize integration-by-parts identities within 
dimensional regularization J8|,|^] . The first thing to realize is that any integral that one can face in computing the on- 
shell fermion propagator can be expressed through eleven basic three loop integrals (topologies) shown in Fig.l. Any 
integral that belongs to a certain topology is considered to be a function of the powers of denominators (both positive 
and negative integer powers are allowed) and one irreducible numerator in each case. For each of the topologies one 
writes down a system of recurrence relations based on integration-by-parts identities that are derived using the fact 
that in dimensional regularization any integral of the total derivative vanishes: 







f — d" 



x propagators) 



(10) 



Here fej is one of the three loop momenta and lj is either one of the loop momenta or the external momenta. Hence 
the starting set of equations consists of twelve recurrence relations for each of the topologies. The fact that not all 
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of these relations are independent is not an obstacle since, when trying to solve them, some relations will vanish 
identically and for this reason will automatically be of no use. 

The next step is to solve the system of these equations. Though there are several ways of thinking about what such 
a solution might be, we prefer to look for the most general one. We then want to construct an algorithm that, for a 
given topology and for any given initial set of powers of propagators, expresses an initial integral through a minimal 
set of "simpler" integrals. The simpler integrals are usually those that either have denominators raised to small powers 
or those that belong to simpler topologies. We then consider these "simpler" , but still non-trivial topologies, write 
down a new set of recurrence relations for them, construct an algorithm that reduces any integral to even simpler 
topologies and continue along these lines until we have an algorithm that completely solves the initial problem in 
terms of a few master integrals. The final set of master integrals is found experimentally. There is no proof that the 
set we find is indeed minimal with respect to integration-by-parts relations in the strict mathematical sense, but for 
practical calculations this set of integrals is sufficient. 




(k) 

FIG. 1. Examples of three-loop quark propagator diagrams corresponding to eleven integration topologies. 

Our solution of the system of recurrence relations shows that it is possible to express any integral which belongs to 
the above topologies through 18 master integrals. Most of these i nteg rals have been calculated in the course of the 
analytical calculation of the electron anomalous magnetic moment []10[ and can be taken from there. It is remarkable 
that a transition from the abelian theory to the non-abelian theory does not result in a significant increase in the 
number of master integrals to be computed, although the number of basic topologies does. As compared to Ref. JlO| , 
we need one additional master integral that corresponds to topology A and we also need one of the master integrals 
of Ref. [[[o| to a higher order in the regularization parameter e. For the QCD wave function renormalization constant 
we also need the constant C\ (see pffl ) which was not computed in pi, because it mysteriously canceled in the 
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calculation of the electron anomalous magnetic moment^. The calculation of these new master integrals is described 
in the next Section. 

There are two principal checks on our solution of the recurrence relations. First, we have computed the three-loop 
anomalous magnetic moment of the electron and confirmed the result of Ref. [^o). Second, the actual calculation of 
the on-shell quark mass renormalization constant Z m has been performed in an arbitrary covariant gauge for the gluon 
field. The explicit cancellation of the gauge parameter in our result for Z m is an important check of the correctness 
of the calculation. There are also some checks related to anomalous dimensions of the heavy quark field in HQET, 
that can be computed combining our result for Z 2 with the known result for Z^ s . We will elaborate on this issue in 
Section V. 

III. NEW MASTER INTEGRALS 

In this Section the calculation of three new master integrals is described. Note that throughout this paper we 
work with the integrals defined in Euclidean space. A pictorial representation of all master integrals is given in 
Fig. 2 in the appendix. The space-time dimension in what follows is parameterized by D = 4 — 2e. We will also use 

oo oo 

C(e) = [7r 2 - £ r(l + e)] 3 , Ck = E l / nk f «r the Riemann ^-function and a 4 = £ l/(2"n 4 ) = Lit (1/2) in what follows. 

n—l n—1 

A. Master integral iio 

Let us consider the master integral Iio. Our basic approach to computing on-shell master integrals is to perform 
a large mass expansion which yields an explicit power series representation for the on-shell integrals. The advantage 
of this approach is that in this way most multi-loop integrations can be performed rather easily, and the resulting 
representation of the on-shell master integral (expanded around the space-time dimension D = 4) takes the form of 
a nested harmonic sum which can be reduced to known mathematical constants. 

This approach has been applied in earlier works |L2| to obtain exact results for several on-shell integrals. It relies 
heavily on an ability to reduce specific classes of nested sums. A similar problem of reducing harmonic sums appears 
in calculations for deep-inelastic lepton-nucleon scattering, where one also deals with an infinite expansion (the light- 
cone expansion instead of the large mass expansion that we employ here) . Collections of summation identities of the 
type that are needed for the present work can be found in the literature Jl^-|T^1 . 

For I 10 we start by taking the squared masses of the two particles inside the loop M 2 to be different from the square 
of the external momenta p 2 — —m 2 . Eventually, we are interested in the result for M 2 = m 2 . 

The expression for Euclidean integral 7 10 reads: 



d g fci d D k 2 d D k 3 
(fci + p) 2 (k 2 + M 2 ) (k 2 +M 2 ) (jfei + k 2 ) 2 (k 3 



J io - / / / 7Z ; „\ 2 ru2 ; *.r2\ n,2 ; ttt\ tz , i \?. it. ; i \i ( n ) 



and we want to perform a systematic expansion in (p 2 /M 2 ). In the present case this can be achieved by making an 
ordinary Taylor expansion in the external momentum p in the integrand 

1 1 VV Y ( 2k 1 -p + p 2 \ l T(a 1 + l ) 

\21cti (b 2 ) a i 2-~S ' I '- 2 I w - ^ -■ ' ^ ' 



[(fci+p) 2 ]"! (fc?)" 1 V k 2 ) r(ai)i 



As a side remark, let us note that for integral I\§ the same large mass expansion can be obtained via a Mellin-Barnes 



representation for the two massive lines |18|. The Mellin-Barnes method can be used successfully here because the 
underlying massless integral with arbitrary complex powers of the massless lines can be evaluated in a closed form in 
terms of a product of Euler T-functions. However the diagrammatic large mass expansion as we use it here (i.e. the 
expansion via a sequence of Taylor expansions in the integrand) can be used also for integrals where the underlying 
massless integral is not known for arbitrary complex powers of the propagators, and the diagrammatic approach to 



3 Moreover, we observed a similar cancellation in our calculation of the three loop slope of the Dirac form factor jnj as well 
as in the calculation of the relation between the MS and the pole quark masses 
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a large mass expansion is therefore more suitable for certain complicated integrals, even though it seems technically 
more involved for simple topologies. 

After performing the Taylor expansion in the integrand in Eq. ( |ll| ) one is left with a 3-loop vacuum (bubble) 
integral with a tensor numerator involving powers of (k\ ■ p). This numerator is simplified after a bubble tensor 
reduction: 



(h-p) n - 

so that Eq.(fL2|) effectively becomes: 

1 1 



T(D/2) 



2"(n/2)! T(D/2 + n/2y for " even ' 
0, for n odd, 



[(k 1+P ) 2 }^ (fc 2 ) 



«i E ( fc? 



p 2 \ 1 £(ai + i) £(ai + i + 1 - D/2)T{D/2) 
fc 2 J r(ai) i! r(ai + 1 - D/2) T(i + D/2) ' 



(13) 



(14) 



The 3-loop vacuum bubble integral can now be evaluated in terms of T-functions 



d^fci d v k 2 d u k. 



{k 2 )^{k 2 + M 2 ) (fc| + M 2 ) [{k x - fc 2 ) 2 ] [(k 3 - fc 2 ) 2 ] 



3 d-8- 2qi r(-3£>/2 + 4 + ai )T*(D/2 - l)r(2 - J/2) 



r 2 (3-£>) 

r(6-2D) 



r(u/2 - 1 + ai)r(Z>/2) 

T{D/2 - 2 + j)r(2 + j - D)T(2 - D/2) 
< (j-l)!r(-3D/2 + 4 + j)r(-l + £>/2) 



(15) 



where ai is any non-negative integer. In this way one obtains an explicit representation for as an infinite power 
series in (p 2 /M 2 ). To obtain the result on the mass shell, where m 2 /M 2 — 1, we have to perform the infinite sum 
over the coefficients of the power series. This task becomes much easier once we expand the obtained representation 
for iio around D = 4. For this purpose the T-functions are expanded as, 



T(n + l + e) =7i!r(l + e) 
where the harmonic sums Sk are defined as 



1 + eSi(n) + — (S 2 (n) - S 2 (n)) + 0(e 3 ) 



n 1 

Sk(n) =J2lk- 



(16) 



(17) 



Expanding in e we obtain the following off-shell expression for 7 10 , written in a form that allows a straightforward 
reduction of the remaining infinite sum as soon as we go on the mass shell 
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The nested harmonic sums S , n ,...,m(*) appearing in this expression are defined recursively via 

Sn ,m(j) 



Sk,n,...,m{i) — ^ ' 
3=1 



r 



(18) 



(19) 



Finally, in the on-shcll limit m 2 — M 2 = 1, the infinite sum over fc in Eq.(|l8|) is performed to yield the following 
result 
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The infinite sums that one needs to go from Eq.(|l8|) to Eq.(|20|) are given in Appendix B 



B. Master integral J5 

Let us now consider a calculation of the constant C\ (see Ref. JlC[| ). The easiest way to extract it is to consider 
the integral I 5 of Ref. ETJ (see Fig.2). We introduce II 1 = fc 2 , II3 = (fci - fc 2 ) 2 , II4 = k\ + 2yk\, II 5 = fc| + 2pfc 2 , II 6 = 
fc| + M 2 , II7 = (fc 2 + h$y + M 2 with p 2 — —m 2 and consider 
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h(m,M)= / nnnnnn (21) 



nin 3 n4n 5 n 6 n7 



for m = M. The result of Ref. 10 reads: 



I 3 1 ( -k z 55 \ 4 4 14 7 2 95 



/ 5 ( m ,H=c( £ ){^ + ^ + ^- y + fj--^--fa3)- r . 2 

+£ (_^_ 44C 3-|^ + l|l + 2 Cl )}, ,22) 

consequently, in order to extract C\ we have to compute I5 to order 0(e). 

For this purpose, it turns out to be useful to derive a suitable integral representation for the vacuum polarization 
subdiagram in I 5 . This representation can be obtained by subsequent application of dispersion relation to "vacuum 
polarization" subgraph and the Mellin-Barnes transformation. Let us first write a dispersion representation for the 
vacuum polarization subgraph: 

d D h 1 7 dA 2 a ./ 4M'\ 1/a - g 

n 6 n 7 7T J A 2 + fcl2i- 2 -r(3/2-e) \ A 2 / ■ 1 ' 

AM 2 

Note that the above expression is written for an Euclidean momenta fc| • We then use the Mellin-Barnes representation 
for the propagator l/(/c| + A 2 ): 



k\ 



TA 2 " = 2^ / + O^D^A 2 )- 1 -^, (24) 



where the integration is performed over the contour C in the complex a plane which goes from complex infinity to 
complex infinity and passes the real axes between the poles of r(— a) and T(l + er). It turns out to be more convenient 
to shift the integration contour to the right, crossing the singularity at a = 0, so that the new integration contour C 
crosses the real axes at e.g. a = 1/2. We obtain: 



A'2 



TA 2 " = A? + hi I + ^(kinx 2 )- 1 -*- (25) 



Let us study two terms in that equation separately. The first term produces a factorized expression where the integral 
over A can be carried out independently of the integration over k\ and £2- The remaining two loop integral is an 
on-shell integral and for this reason can be easily done. The result of this calculation reads (we take the limit m = M): 

r t s „t \ f 1 3 1 / 7T 2 55\ 8 A 2 43 

In order to compute the contribution of the second term in Eq.(^5|) we insert it into dispersion integral Eq.(|23|) and 
integrate over A. We obtain: 



h b (m, M) = i~ 2l XZ-e) I d ^~^ 1 + ^ {* + £ ' f ~ £ ) ( 4M2 ) 
d^fcid^fca 



2\-E-CT 



n!n3n4n 5 



(n 2 ) CT , (27) 



where II2 = k%. 

To proceed further, we need to perform the integrals over k\ and k^. The easiest way to do that is to utilize the 
integration- by-parts identities. It turns out that there is an identity that allows to remove either II3 or II4 or II5 in 
which case the integrals become computable in a closed form for arbitrary a. To write down this identity, we introduce 
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. d D k!d D k 2 , , 

I(ar,a 2 ,a 3 , ai ,a 5 ) = / f ^ f ^ f ^. (28) 



The general form of the identity is 19 



I{{ai})= ? 1 



1 f 2(L> - 1 - a 3 - ai)(D - 1 - a 3 - a 2 ) 



2(2£> - 2 £° =1 a* + a 4 + a 5 ) I (£>-l-2a 3 ) 
(£>-l-a 3 -a 2 ) i+5 _ i _ _ (fl-l-a 3 -ai) + 
(D - 1 - 2a 3 ) (£> - 1 - 2a 3 ) 

+a 3 3+(4- - 5~)(2- - I") - (a 2 - 1)5" - (o x - 1)4~ }/({a<}), (29) 

where as usual 1 /(fli, ■■■) = I(a\ ± 1, ■■■) and similar for other operators. 

It is now easy to see that if we apply that recurrence relation to the integral 1(1, —a, 1, 1, 1), we immediately obtain 
a set of integrals that can be easily computed since the recurrence relation removes either one of the massive lines 
(in which case there appears a massless two point subgraph) or the line II 3 (in which case the integral becomes a 
product of two one loop integrals). Therefore, we end up with the expression for 1(1, —a, 1, 1, 1) written in terms of 
T- functions: 



l( I . -n. 1 . I . I ) = - ^^S 2 (1, -a, 1, 1) - is 2 (l, 2, -1 - a, 1) 



1 

2 + 2(7 - 4e 

is 2 (2, -1 - a, 1, 1) - is 2 (2, -a, 0, 1) + ( (l - e + ±a) S 2 (l, -a, 0, 2) 

\-e) S 2 (l, 1,-1 - cr,2) + (-2 + 2e<7 + 6e - 4e 2 - cr) &(!, 1) 5i(-<r, 1))] , (30) 



where the functions S 2 and Si are defined as follows: 



S 2 (ai,a 2 , a 3 , a 4 ) = O x (ai, a 2 ) Si(a x + a 2 + a 3 - 2 + e, a 4 ), 
„ , , T(2 - £ - ai)r(2 - e - a 2 )r(tn + a 2 - 2 + e) 

Ci(ai,0 2 J — =7 rp-; rpr-j r , 

T(ai)r(a 2 )r(4 - 2e - a\ - a 2 ) 
a , v r( ai +a 2 -2 + e)r(4-2e-2ai -a 2 ) 

Si(ai,a 2 ) = , - , — r . (31) 

1 (a 2 )l (4 — 2e — a\ — a 2 ) 

We use these expressions in Eq. ( ^j ) and integrate over a by taking the residues located to the right of the integration 
contour, i.e. at the points a = n, n + e, n + 2e, where n > 1 is an integer^. Finally we arrive at the representation 
for 1$ written as an infinite series: 

oo / 2 \ n 

I 5b (m,M) = J2dn(^) ■ (32) 

The expression for the coefficient d n is too long to be presented here; the essential point, however, is that d n can 
be expressed through harmonic sums and the summation for m = M can be carried out in a way similar to the one 
described in the previous subsection. Our final result for I§b (m,m) reads: 

h b (m, m) = -^tt 4 - 2C(3) - ^tt 2 + 26 
45 3 

+e U + f + 4^ 2 log 2 - | ^ - 26C 3 - ^ - ^fc) • (33) 



Combining the results for I$ a and /5b and comparing with Eq.(|22|), we obtain the constant C%: 

C x = Cs - - — 7T 4 + — tt 2 C 3 + 2^ 2 log 2 - — C 5 . (34) 
SJ 3 180 12 ^ 6 4 S5 v y 



4 It is interesting to note that different poles in a correspond to contributions of different expansion regimes if I5 is expanded 
in m/M according to the rules of the large mass expansion (for a review see e.g. pp|). 
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C. Master integral lis 



For integral lis we proceed as for Iiq and perform an expansion in p 2 /Ad 2 of a more general Euclidean integral 



'18 



d^fci d D k 2 d D k 3 



(k x +p) 2 (fc 3 + P) 2 (h - k 2 ) 2 (k 3 - k 2 ) 2 (k 2 + M 2 ) (k 2 + M 2 ) (k 2 + M 2 ) 



(35) 



As is the case for I\q , the expansion of lis can be performeoj^] by making an ordinary Taylor expansion in the external 
momentum p in the integrand of Eq. (|35"|), i.e. by expanding both l/(fci +p) 2 and l/(k 3 + p) 2 in small p as in Eq.(fl2"|). 

The resulting integral is formally a 3- loop vacuum bubble integral that, after the necessary tensor reductions and 
partial fractioning of massive and massless lines, can be evaluated in terms of T- functions using Eqs.(|67|-[7ll) in the 
Appendix A. In this way one obtains an explicit representation for Ii$ as an infinite series in p 2 /M 2 with coefficients 
being a nested sum over a product of T- functions. After expanding the expressions in s as in Eq. (|l6|) and writing the 
nested sums in a form that allows a straightforward reduction on the mass shell we obtain a fairly simple expression 



/ 18 = 



C(e) 

M 2+6 S 



^ ( M 2 



k=l 



9 C3 _ „ (2 
k 2 k S 



4C 



si(k) sm 

k 2 fc 4 fc 3 



| ^i(fc)^(fc) ^(k) 



k 2 



k 2 



0(e) 



(36) 



The on-shell result is now easily recovered by putting m 2 
Eqs. @p|). We finally obtain 



M 2 



1 and using the expressions for the infinite sums 



hs = C(e) 



2tt 2 ( 3 



5Cs + 0(e) 



(37) 



Having derived the last unknown contribution for the three loop on-shell two-point master integrals, we are able 
to write down a complete set of master integrals for the three loop on-shell two point functions without unknown 
constants. The list updates the result of Ref. JTol and can be found in Appendix C. 



IV. THREE LOOP QCD RENORMALIZATION CONSTANTS IN THE ON-SHELL SCHEME 

Performing an explicit computation along the lines outlined in the previous Sections, we are able to derive an 
expression for the QCD on-shell renormalization constants up to the three loop order. We work here in a general 
covariant gauge with a gauge fixing parameter £; the gluon propagator (omitting the color indices) therefore reads: 

D^P) = -^U»-^). (38) 

In our previous article || we have given a relation between the pole and the MS mass, from where the value of Z m 
can be determined. Here we give the result for that renormalization constant explicitly. We write: 

Z m = 1 + a Q C F (- 1 - j^^j + a 2 C F z£> + a 3 C F Z^ , (39) 

where 

_ ai Q) r(l + e)m- 2e 



5 We should note that a large mass expansion of integrals 7ig and Jio is rather special since no contributions arise other than 
the Taylor expansion of the whole integrand in p. However, for other types of integrals without a through-going massive line 
the combinatorics of the large mass expansion gives rise to several more subgraphs that must be Taylor expanded. 
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and af is the bare coupling constant and Z$ and Z„' are the two and the three loop coefficients respectively: 
Z%> = C F df + C A df + T R N L df + T R N H df , 

Z%> = G% df + C F C A df + C\ df + C F T R N L df + C F T R N H df + C A T R N L df 

+C A T R N H df + T 2 R N L N H df + T 2 N 2 H df + T R N 2 df . (40) 

Here Cf and Ca are the Casimir operators of the fundamental and the adjoint representation of the color gauge 
group (the group is SU(3) for QCD) and T R is the trace normalization of the fundamental representation. Nl is the 
number of massless quark flavors and Nh is the number of quark flavors with a pole mass equal to to. 
Our results for the coefficients d^ are: 

,(2) 9 45 199 3 , 1 2l 5 2 
^ = 32^ + 64^ + 128-4 C3+ 2 7r ^ 2 ~ 16^ 

( 677 33 > „ „ 9l n 55 9 7 A 1, 4 „ 
+£ ( 256 - T Cs + ^ l0g 2 - ^ l0g 2 - 32 * + 40 " " 2 l0g 2 " ^ 
,(2) 11 91 605 3 , 1 2l „ 1 2 

/ 3799 13 3 2 , „ 1 a, 2„ 19 2 7 4 1 4 „ „ 
+£ (" W + T C3 - 2^ l0g 2 + r l0g 2 + 48* " 80* + 4 l0g 2 + 6 ° 4 

,( 2 ) 1 7 45 1 2 ^279 , 7 2 
d3 + 32 + 12^ +e UT + C3+ 24^ 



12^ 


+ e 


V"6T 






'463 


6 


-( 






457 


9 

+ 16 


K- 


512 



,(2) 1 7 69 1 2 / 463 7, 2l „ 5 2 
^^8? + 16^ + 32-6^ +£ U-2 C3 + 7rl ° g2 -6 7r 

(3 ) 9 63 1 / 457 9 , 3 „ 15 A 1.4225 

"i 



128e 3 256e 2 s \ 512 16 s 8 b 64 j 3072 

-> 2 + & + s & + '- 2 loe2 + r 3 *" 2 - si- 2 - 1^ 4 - 1 b » 4 2 - 3 - 

, (3 ) 33 1039 1 /2563 53 53 2l 61 2 

^ =l^ + 768? + H^12--32 C3 + 48 7r ^'-ge" 

144929 19 , 2 2459, 45, 223 2l „ 221 2l 2 „ 

H C37I" C3 H C5 + 71" log 2 7T 2 log 2 2 

9216 16 s 96 y 16 s 36 8 72 B 

1477 2 4639 4 167, 4 „ 167 

7T 2 H 7T 4 bg 4 2 a 4 , 

576 8640 144 s 6 ' 

l(3 ) 121 167 1 / 81797 11 11 2l 11 , 

* = -5^-T08? + irT0368 + T6 C3 " 24 7r l0g2+ 72 * 

2188703 51 2 3059, 65, 313 2l „ 

1 C37T H C3 C5 7T 2 log 2 

62208 64 s 288 s 32 s 72 6 

11 2 , 2 „ 553 2 3667 4 11, 4 „ 44 

H 7T 2 log 2 2 H 7T 2 7T 4 H log 4 2 H a 4 , 

9 B 3456 17280 18 8 3 



,(3) 3 77 1 / 415 1 1 21 „ 7 2\ 2911 169, 

d ' --^-W^ + -e\-WA + A^-r l0g2+ 48" J - 2304 + ^4 C 
29 2l „ 8 2l 475 2 371 4 4, 4 „ 32 

7T 2 log 2 + -7T 2 log 2 2 H 7T 2 7T 4 + - log 4 2 H 4 , 

9 B 9 8 288 2160 9 B 3 
,(3) 3 77 1 / 631 1 , 1 ,, 1 A 8743 71 , 

4 = "^- 192? + 7 ("384 + 46- 3* l0g2+ 3" J~ 2304 + 12 C3 
67 2 , n 5 2, 2„ 425 2 161 4 4, 4n 32 
-36^ l0g 2 + 9* l0g 2 + 432^ ^ 2160* + 9 l0g 2 + i" 4 ' 
, (3 ) 11 877 1 /3125 1 2l „ 7 2 \ 317381 41 , 
d « = 72^ + 864^ + UW + 6" l0g2+ 72" J + T555T + 144 C3 
29 9 , „ 4 9l 2 „ 73 2 29 4 2 , 4 „ 16 
+ 18* l0g 2 - 9* l0g 2 + 108* + 432* - 9 l0g 2 - T fl4 ' 
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,( 3 ) 11 877 1 /502 1 2l „ 13 ,\ 473549 1 2 1345, 
4 = 72^ + 864^ + H^ + ^ l0g2 -36^ J + T5552- + 8^ ~ UA^ 
5, H5 2 , n 5 2l 2 „ 707 2 1 4 2, 4n 16 

( 5 H 7T log 2 7T 2 log 2 2 7T 2 + — 7T 4 - - Og 4 2 fl 4 , 

8 s5 18 8 18 8 144 54 9 8 3 ' 
,(3) 1 17 1/152 1 2 A 2032 1 7 2 13 2 

4 =-— -—, +^-^ + 18- J " M3 + " 3" lo S 2+ 27"' 

1 / 385 1 2 \ 5441 53 A 2 2 , „ 79 2 
d 9 =-053- 17553 + 11-^4 + 9' J" gra + IS 6 " 3" l0g2+ 135"' 



18e 3 


54e 2 


1 


17 


36e 3 


108e 2 


1 
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4o } = -^3 ~ 777^2 + 7 ( -^77 - - ^ ) - — - -C 3 - -^r 2 . (41) 



1 / 223 1 2 \ 2687 19 17^ 
36e 3 108e 2 ' e V 324 ~ 18 ^ J ~ ~972 ~ W 3 ~ 54 ' 

oo oo 

We have used here (fc = l/n k for the Ricmann (-function and also a 4 = ^ l/(2"n 4 ). 

n— 1 n— 1 

It is more cumbersome to present the result for the wave function renormalization Z 2 since, as we mentioned already, 
it becomes gauge dependent at O(af). We parameterize Z 2 as: 

Z 2 = 1 + a C F (-^ - y~2^) + a o°F Z 2 2) + a l C F Z 2 3 \ (42) 

where 

4 2) = Cf A (2) + C A / 2 (2) + T R N L + T R N H /f \ 

Z (3) = C 2 A (3) + C F C A / 2 (3) + Cft / 3 (3) + C F T R N L /| 3) + C F T R N H / 5 (3) + C A T R N L / 6 (3) 

+C A T R N H + T 2 R N L N H f (3) + T 2 R N 2 H / 9 (3) + T 2 N 2 f$ . (43) 

Our result for the coefficients /j^ reads: 

,(2) 9 51 433 3 2, o 13 2 

-3^ + 647 + 128-2 C3 + 7r l0g2 -16^ 

/211 147, 23 „ „ ,„ 89 , 7 4 , 4 „ „ \ 

+£ ( 256 - — C3 + T* l0g 2 - 2 * bg 2 - 32 * + 20 * " l0g 2 ~ 24 ° 4 j ■ 
r(2) 11 101 803 3, 1 2l „ 5 2 
/2 --3^-647-128 + 4 C3 -r l0g2+ 16- 



4241 129 , 23 2l „ 2 „ 41 , 7 4 1, 4 „ „ \ 

1 C3 7T og2 + 7T 2 bg 2 2 H 7T 2 7T 4 + - log 4 2 + 12a 4 

256 16 ^ 8 8 8 48 40 2 8 *) 



( 2 ) 1 9 59 1 2 / 369 . 3 2 



f * =8? + lfe + 32 + T2 7r +£ U4 +C3+ 8 7r 



(2) 1 19 1139 1 2 / 20275 , „ 2l „ 19 2 
4 -4? + 48l + W-r +£ (l72y- 7C3 + 2 ^ 1 ° g2 -l2 7r 



,(3i 9 81 1 / 1039 9, 3 2l „ 39 2 \ 10823 1 , 

' fl = -^- 256? + H-M2- + 8 C3 -4 7r 1 ° g2 + 64 7r J-^072- + 8 C3 " 

187 A 5 , 685 ?1 „ „ o, 2 „ 7199 2 41 4 5 , 4 „ _ 

C 3 + log 2 + 3tt 2 bg 2 2 7T 2 7T 4 bg 4 2 - 10a 4 , 

32 S3 16 S5 48 8 8 1152 120 12 8 4 ' 

,( 3) 33 1217 1 / 14887 27, 53 2 , „ 331 2 
/2 = i^ + 768? + e( v T536--Y C3 + 24 7r 1 ° g2 - 192 7r 
150871 45 2 9941 145 A 2281 2l 
+ ^216- " 16 C37F - T92" C3 + T6- C5 + -m* l0g2 
499 2 , 2 U69 2 20053 4 319 , 4 319 

bg 2 - + 17280^ ' 144 bg 2 " -° 4 ' 

,( 3) 121 1501 1 / 55945 173 , 11 2l „ 55 2 



576£ 3 864e 2 e \ 5184 128 12 b 96 1080 
3 7T 4 1 \\ 2551697 127 2 14371 37 

" s ' ! ~ 256 C3 + 4320 ~ 768 J J ~ 62208 + ~T2^' K + 576 ^ ~ Y^ 5 
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271 2l 391 2l 2 1073 2 10811 4 349 , 4 349 

tt 2 log 2 H 7r 2 log 2 H tt 2 7T 4 H log 2 H a 4 

36 s 144 s 3456 23040 288 6 12 

. / 13 1 2 13 7 . 1 2 17 4 

+ H"768 + 144 C3 ^ " 256 C3 + 384 Cs ~ 256^ + 27648^ 
r(3) 3 103 1 / 351 3, 2 2i 23 2 \ 3773 413, 

/4 = -32^-192^ + H _ T28 + 4 C3 "3 7r bg2 + 48* ) " 2304 + H Cs 
58 2, „ 16 2l 2„ 905 2 733 4 8, 4 „ 64 
"IT bg 2 + <T l0g 2 + 288^ 2 " 2160* + 9 bg 2 + T ^ 
,(3) 3 91 1 / 1525 , 2 2l „ 19 2 \ 78967 5273 , 
h ^~T6T^T92^ + -e{-^M +C3 -3 7r bg2 + 24 * ) 6912 + W Cs 
SI 2, „ 5 2l 2„ 865 2 1 3 7 4 7, 4 „ „„ 
"IT bg2 + 6^ bg 2 + 648* " 720* + 6 bg 2 + ^ 
„ f31 11 1069 1 /550 1 1 2l „ 1 2 \ 416405 1 4 5 , 29 2l 

^ 1 -ra^ + sei^ + nir-^ + r 1082 -i8^J + 15552 - le^ + T' log 

8 2 , 2„ 7 9 29 4 4, 4 „ 32 
— 7T 2 log 2 2 H 7T 2 + 7T 4 - - bg 4 2 a 4 , 

9 B 18 216 9 6 3 ' 

,( 3 ) 1/15 , 1 \ 1 /353 , 1 \ 1 /503 1 2 , „ 59 o 1, 35 

^ - ?U-^J + ?( v 288 + ^64j + n^ + 3 " ^72* " 2°» " f 576 
32257 11, 2 13301 , 15, 521 „ 1 2l 2 „ 13567 , 
+ -64T + 48 Cs * - ^76- C3 ~ 16 Cb + 36 ^ bg 2 " 3 bg 2 " W* 

+ ^ 4 -^log 4 2-16a 4 + ef^^C3 
72 3 & ^V1728 24 s3 

,.(31 1 7 1 / 31 1 2 \ 1168 4 2i „ 7 2 

'» - -12? ~ i8^ + e hr + D * J - + 4C3 - r log2 + 6* > 

/f = * * +±(-^ + -« 2 ) - — + 7Ca - ^ 2 log2 + — 7T 2 , 
J9 12e 3 36e 2 e \ 54 9 / 648 w 3 8 10 ' 

.(3) 1 23 1 / 325 1 2 \ 4025 19 . 23 2 

/in = 5" o + 7T C3 1" ■ (44) 

10 36e 3 108e 2 e V 324 18 / 972 18 s 54 v ' 



Here the gauge parameter £ is defined in Eq. ( |38| ) and we again stress that the three loop contribution to Z 2 is gauge 
dependent. This fact is not too surprising by itself, since the wave function renormalization constants are known to 
be gauge dependent in general (as one can see for instance in the MS scheme). The on-shell renormalization constant 
computed in dimensional regularization, however, unexpectedly turned out to be gauge independent in two first 
orders in perturbation theory. This fact, observed for the first time in |Q, resulted in "all orders gauge independence" 
conjecture by the authors of that reference. Our calculation shows that things "return back to normal" in the 
third order of perturbative expansion, where gauge parameter dependence explicitly appears. Let us note that the 
"strongest" gauge dependence appears in CaTrNh color structures. We note in this respect that the JVjj-dependent 
divergences permit a strong check related to decoupling of heavy quark loops from HQET. We discuss this issue in 
detail in the next Section. 

V. HQET WAVE FUNCTION RENORMALIZATION 

There is an interesting application of the results derived above which also serves as an additional check on its 
correctness. Namely, the on-shell wave function renormalization constant and the known MS wave function renor- 
malization constant can be combined to derive the heavy quark wave function renormalization in HQET. To two 
loops that has been done in Ref. 0. Moreover, the HQET wave function renormalization constant^ can be shown to 



6 When we refer to HQET wave function renormalization constant we mean the MS wave function renormalization constant 
in HQET. 
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exponentiate , i.e. it can be written asQ 

^HQET = e (^C F x 1 + (^) 2 C F (CA^2+TRAr i X3) + (^) 3 C F (cix4+CAT R Af I/ ^5 + A'£T|x 6 +C F Af L T ii:C7 )+0( C[ f)) ^.v 

For this reason, only few genuinely new color structures appear at the three loop level and the Cp and C f Ca 
coefficients in Z^ ET are completely determined by the two loop expression for Z^ ET . This remarkable feature 
provides an additional strong cross check on our result. 

Let us first explain why is it at all possible to determine the HQET wave function renormalization from the MS 
and the on-shell wave function renormalization constants in full QCD. The explanation is apparent if one asks what 
HQET implies for the computation of Feynman diagrams. In fact, all underlying assumptions of HQET are satisfied 
if we try to compute the quark self energy in the vicinity of the particle mass shell. If we denote the off-shellness by 
—A 2 = p 2 — m 2 , the HQET Lagrangian gives a prescription of how non-analytic contributions in A to the quark self 
energy should be computed. Such a computation should clearly be supplemented by a calculation of the contribution 
coming from momenta region k A, which is analytic in A. This contribution is provided by the Taylor expansion of 
the quark self energy in A in full QCD. At this point the computation is reduced to the computation of the on-shell 
integrals in full QCD and this is precisely the piece of work we do when we compute the on-shell wave function 
renormalization constant in full QCD ( i.e. we compute the quark self-energy operator to the first non-trivial order in 
A). The on-shell wave function renormalization constant in full QCD is both ultraviolet and infra-red divergent. The 
ultraviolet divergences can be removed by means of the MS renormalization; the infra-red divergences then cancel out 
once the HQET contribution to the quark self energy is computed since the off-shell self energy is infra-red finite. We 
therefore see, that once the infra-red divergences of the on-shell Z% are isolated, they should match the ultraviolet 
divergences in HQET and this is the way how the HQET wave function renormalization constant and the HQET 
anomalous dimension of the heavy quark field are determined. 

We now turn to the computation of this HQET quantity. In order to isolate the infra-red divergences, we consider 
the ratio of the MS and the on-shell wave function renormalization constants in full QCD. Before giving the complete 
result, let us first discuss what we expect to happen to the contribution of heavy fermion loops. 

To begin with, it is well known that the contribution of heavy quark loops in HQET vanishes identically due to the 
fact that the anti-particle pole is removed from the heavy quark propagator. Therefore, a quark- antiquark pair can 
not be created and this is the reason why the heavy fermion loops are disregarded in HQET. Hence, it is impossible 
to obtain the color structures proportional to Nh out of any HQET computation. As it follows from the preceding 
discussion, this should imply that in the ratio Z^ s /Z® 5 , computed in full QCD, any A#-dependent structure should 
become finite since otherwise it would be impossible to remove the remaining divergence by performing a HQET 
calculation. Note, that this requirement on A^y-dependent structures provides a check on the gauge dependence of 
Z 2 computed in this paper since the color structure CaTrNh is gauge dependent. By explicit computation of the 

ratio Z^ s /Z® s we find, however, that a divergent pieces proportional to Nh remain in this ratio. One should recall 
at this point that the decoupling of heavy degrees of freedom in field theories within the minimal subtraction scheme 
is not automatic. In order to observe the decoupling of heavy particles, one should express all MS renormalized 
parameters of the theory in terms of their effective parameters whose renormalization scale dependence is governed 
by only the light degrees of freedom. In QCD this is done using well established decoupling relations 0j22|. Since 

the ratio Z^ s /Z° s depends on the coupling constant and the gauge parameter, one should take into account the 
decoupling relations for these two quantities in order to obtain a proper low energy result. Once ots H+N and 
£(n h +n l ) are ex p ressec i m terms of ai NL ^ and £ ( Nl \ one observes explicitly that all divergences in the A#-depending 
terms disappear and one is left with the expression that can match the structure of divergences in HQET. Let us 
emphasize that this argument alone proves that the on-shell wave function renormalization constant Z2 in full QCD, 
that we compute in this paper, should be gauge dependent ! 

After this discussion, we are ready to present our result for the HQET anomalous dimension of the heavy quark 
field in the MS scheme. The expression for the wave function renormalization is then easily derived from the condition 
(Zp/Z 2 os )Z 2 HQET = finite. We define: 



dlog^ 2 HQET 
d log jj? 



7HQET- (46) 



7 We are not aware of any systematic discussion of this point within HQET. The clue that there is such an exponentiation, can 
be taken from earlier studies of the eikonal approximation in QCD [E3[. We thank A.G. Grozin for discussions on this point. 



14 



Our result for the HQET anomalous dimension 7hqet then reads: 



where 



7i = 


c F 




72 = 


c F 


C A ( 


73 = 


c F 


c\{- 


-+ 


H 


69 




2048 



-c a t r n l 



7HQET 



1 

2 

19 

"24 
3 

16 
3 
512 
/1105 



7T 



72 I — 



73 (_±) +0 (a*), 



32 s 64 s 



C3 







19495 
27648 " 
5 



1024 



1 

360^ 

■3 



TlNl 



V6912 + 4^ 3 + 256^ 



379 
2048 
5 

108 



15 „ 1 

(3 H 7 

256 1440 



CfT r N l 



51 
64 



(47) 



(48) 



Here a s is the MS coupling constant in the theory with Nl massless flavors and the gauge parameter £ is defined 
in Eq.([38|). The above set of equations gives the result for the MS anomalous dimension of the heavy quark field in 
HQET up to three loops. The two-loop result was obtained directly within HQET in ^4j. The iVi,-dependent terms 
in Eq.(|4"8|) agree with preliminary results of the ongoing three loop calculation within HQET.^J 



VI. SUMMARY OF QED RENORMALIZATION CONSTANTS 

In this Section we summarize the results necessary for the on-shell renormalization of QED up to the three loop order 
and also demonstrate that Z2 in QED is gauge parameter independent in a general covariant £ gauge, if dimensional 
regularization is used for both ultraviolet and infra-red divergences. This result for Z2 was previously mentioned in 
; in other regularization schemes the gauge dependence of Z 2 was studied in . While we do not add anything 
new to this question in essence, we think that the derivation presented below is relatively simple and transparent and 
so we decided to include it into our discussion for completeness. 

Let us start with deriving an equation which determines a dependence of the fermion propagator on the gauge 
parameter £ in a general covariant gauge. Such an equation is easiest to derive in the path integral formulation of 
QED. In a general £ gauge the expression for the two-point electron Green function reads: 

r — - / dxB(x) 2 r 

(T$(*M0)> 4 = A^ 1 / [dB] e 2 VJ I [d$\ [#] [dA]S(d^ - B)j(x)^(0) e iS , (49) 

where rj = 1 — £, S is the standard QED action 

S = Jd+x (-^F^ +^(id-m^- eofa^A^j , (50) 

and N is the normalization factor: 

r — - / dxB(x) 2 r 
N= [dB] e 2 V J / [d^][d^][dA]*(fl p Ap - B) e iS . (51) 

If the integration over B is performed first, one recovers the familiar expression for the path integral representation 
of the fermion two-point function in a covariant gauge. We now perform a gauge transformation in the integral over 
matter fields: 



8 We were informed by A.Grozin that he is calculating the three-loop HQET anomalous dimension using a recently constructed 
algorithm for 3-loop calculations in HQET J5f|. He has presently obtained all 3-loop terms except for the CfC\ one. The 
result coincides with Eq. (Eg). 
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A„^A tt + d„f, if> -» e~ ie °^, (52) 

-)2 



and choose f(x) to satisfy the equation d f = B(x). We obtain: 

dxB{x) 2 



(T^(x)^(0)) 4 = (T^MO))^ 



-i 



dB] e 2f l 



J[dB(x)} exp ^ J {dx)B 2 + te J 



(dy){G{x-y)-G{-y))B{y)\ (53) 



where G(x) is the solution of the equation d 2 G(x) — 5^(x) and (TyS(x)^(0)).B=o is the fermion Green function in 
the Lorentz gauge d^A^ = 0. 

We can now integrate over B{x) and a convenient way to do this is to perform a Fourier transform and integrate 
over its Fourier components B{k). We then obtain: 



(dk)J k (x)J- k {x) 



(Tty(a;)V>(0)) £ = e 2 J ■ ■ <TV(aM0)) B =o, (54) 

where (dk) — d D k / (2tt) d and </& is defined as: 

■M*) = ^T (e^-l). ( 55 ) 
By differentiating with respect to £ in Eq.(|5^), we obtain 

^(T^(x)^(0)h = -*■§/ (dk)Mx)J- k (x) (T^(x)^(Q)) c . (56) 
In momentum space this equation reads: 



( d l) a _,_ , 2 /" (<**) 



— 5 F (p) 5 = ie y — S F (p - q)t - ie j -j^Spip)^ (57) 

We now analyze Eq.(|57j) close to the mass shell where we approximate the fermion propagator by: 

Sf(p) « -fi-. (58) 
p — m 

The fact that the right hand side of Eq.(|57|) does not develop a double pole in the limit p 2 — > m 2 implies that dm/d£ 
is zero, in other words the pole mass of the fermion in QED is gauge independent. If we equate the single pole 
contributions on both sides of Eq.(|57|), we obtain: 

±lo g Z 2 = -ieljW. (59) 

Note that this result is derived under the assumption that the integral 

(dq) 



S F (p-q) t (60) 



does not develop a pole if p 2 —> m 2 and one can make sure by power counting that this is a non-trivial statement. 
The correct approach, which matches perturbative calculations, is to take the limit p 2 — > m 2 first assuming that the 
integrals remain finite because of the regularization. In this sense the integral in Eq.(|60|) does not develop a single 
pole and can be neglected. We now return to Eq.(|59|) and note that in dimensional regularization the integral over k 
is zero. Hence we obtain: 

-|logZ 2 = 0. (61) 
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We thus proved that the on-shell fermion wave function renormalization constant in QED is gauge parameter inde- 
pendent. 

Let us now turn to the presentation of the three-loop QED renormalization constants. We consider QED of a single 
electron; that implies that Nh = 1 and Nl = should be used in general formulas of Section IV. 
We write the mass renornialization constant Z m as 

(i) 

where cxq is the bare QED coupling constant. For the expansion coefficients Zm we obtain: 

z.w- -— L_ 



4e l-2e 

m 13 73 475 3, 1 . 23 , 

Z » = 3^ + 6te + 128 " I 6 + 2* l0 ^--^ 

/2529 47 9 9 o 245 9 7 A 1 A 

^ ( W - T Ca + ^ log 2 - . 2 log 2 2 - -. 2 + -. 4 - - log 4 2 - 12a 4 

„ (% s 221 5561 1 / 154445 13 , 17 2 , 391 2 

Z « = - ^ + - 41472 + T6 Cs " 24* bg2 + 576* 

5, 89 5l „ 65 2l 2 „ 
Cs + ^5 + ^ log 2+-^ log 2 2 

23 

log 4 2 + — a 4 . (63) 



221 


5561 




1152e 3 


6912e 2 




3489365 


979 
4320 


719 


248832 

5783 2 

7T - 

17280 


' "72~ 
72 



3 

The wave function renormalization constant is parameterized as 



^ = l + Ey^(l + e)J Zjp, (64) 



where 



3 


1 




~4e 


l-2e' 




17 


229 


8453 


32e 2 


+ 192s + 


1152 


■(? 


7T 2 log 2 - 


419 . 

7T 

96 



Z\ 2> = H 1 7T 2 - -C 3 + 7T 2 log 2 

- qi,-2 mo- nco /^g 2 

86797 203, „ 2l 2 „ , 4 „ 7 
C 3 - 2tt 2 log 2 2 - log 4 2 - 24a 4 + 



6912 8 ° ° 20 

^ (3 ) 131 2141 1/17, 17 2 , „ 935 , 116489 \ 

Z\>= ? r + - — C 3 7T 2 log2 + 7T 2 

2 384e 3 2304e 2 e \ 8 12 5 576 13824 J 

2121361 1 „ 9 2803, 5 A 1367 „ 
"^294T + 8 C37r + 144" C3 " 16 C5 + 144* bg2 

23 2l 2 „ 197731 9 383 4 3, 4 „ ,„ 
+ lT l0g " 2 - 1*840 * ~ 720^ + 4 bg 2 + 18a4 - ^ 

We explicitly see that the QED result for both the mass and the wave function renormalization constants does not 
depend on the gauge parameter, in agreement with Johnson-Zumino identity Eq.(|6l|) and in variance with the QCD 
calculation of the previous section. 

Finally, for completeness, we give the relation between the bare ao and the physical QED coupling constant a 
defined in the on-shell scheme (i.e. through the photon propagator at zero momentum transfer). The result reads [Q: 

^ = 1 ~H J' (2- £ )(l- 2g )(l + 2e) (1 + £(7 ' 4£)) ( v ) ' (66) 

and completes the set of the renormalization constants required for the renormalization of the QED up to three loops 
in the on-shell scheme. 
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VII. CONCLUSIONS 



We have presented a calculation of the three loop on-shell renormalization constants both in QED and QCD. Explicit 
result for the mass and the fermion wave function renormalization constants are derived. Dimensional rcgularization 
is used to regulate both ultraviolet and infrared divergences. 

Our calculation represents an important step towards three loop calculations with heavy on-shell quarks in QCD 
(semileptonic 6 decays, NRQCD, etc.). Furthermore, these results are important for applications in QED where the 
on-shell scheme is clearly superior as compared to any other scheme. 

The use of dimensional regularization has some interesting consequences for the gauge dependence of the fermion 
wave function renormalization constant. First, in QED turns out to be gauge invariant to all orders in a. In 
QCD, only the first two orders of the perturbative expansion of Z-i in a s are gauge independent; unfortunately, the 
dependence on the gauge parameter explicitly appears in the third order of the perturbative expansion. 

As a byproduct of this analysis, we derived the anomalous dimension for the heavy quark field in HQET. 
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APPENDIX A: ELEMENTARY MULTILOOP INTEGRALS 



In this appendix expressions are given for a number of integrals in dimensional regularization for which exact 
analytic results are known. We present these below in Minkowski space. 
For a 1-loop massive bubble integral one has Q 

d D p 



{p 2 - m 2 + ie) Ql (p 2 + ie) Q2 



i ^/Vl)-— (mT /2 -— ^ + °T D J?3?J 2 " ■ (67) 



A 1-loop massless propagator-type integral has the simple form 

d D p 



! 



(p 2 + ie) Q i [(p + Q) 2 + ie] a * 



i 1 - -k d ' 2 ( Q 2 )D/ 2 -^-^ r ^ 2 ~ ~ ^)rK + a 2 - g/2) _ 

r(ai)r(a 2 )r(£) - ai - a 2 ) 

A compact expression for this integral with a general tensor numerator can be found in Ref . j|] . For a 1-loop on-shell 
propagator-type integral one derives 

d D p 



(p 2 + ze) Q i {p 2 + 2p ■ Q + ie) a i 

i n n/2 {Q 2 )D/ 2- ai - aH _ irai - a2 r( ai + a 2 - D/mP - 2a, - q 2 ) 

r(a 2 )r(I? - at - a>2) 

and for a two-loop bubble integral with one massless and two massive lines one finds (26) 

d D p d D k 



(p 2 - m 2 + ie) Q i [{p + k) 2 + ie] a ^ {k 2 - m 2 + ie) a z 

Tr D (m 2 ) D ~ ai ~ a2 ~ a3 ( i) 1 -"!-^-^ + a i + Q 2 + Qft) 

r(ai)r(a 3 ) 

Tj-D/2 + ai+ a 2 )T{-D/2 + a 2 + - a 2 ) 

X r(D/2)r(ai + 2a 2 + a 3 - £>) 



(70) 



Also for this integral with a general tensor numerator a compact expression is known [f27f . 

Several more simple cases follow by using these expressions recursively as the powers ai,a 2 and a 3 , are allowed to 
be non-integer, possibly containing D. In this way one can obtain, for instance, 

d D p d D k 



(p 2 + ie) ai (k 2 - m 2 + ie) a ^ (k 2 + ie)" 4 [(p + A:) 2 + ze]" 3 

7r D (m 2 ) D " 0l " Q2_Q ^ a4 ( i)l-°»i-«a-«3-a4 r ( ai + Q 2 + ^3 + «4 - ^P) 

r(ai)r(a 2 )r(a 3 ) 

r(ai + a 3 - D/2)Y{D/2 - a x )r(£)/2 - a 3 )r(D - a x - a 3 - a 4 ) 



r(£»/2)r(D - ai - a 3 ) 



(71) 



APPENDIX B: SOME HARMONIC SUMS 

In this Appendix we present a partial list of harmonic sums which we used in this paper. These results can be 
found e.g. in m4 



19 



oo 1 

^- = 5 fc (oo) = a (*>i) 



E 

2 = 1 
OO 

E 



Si(i 



Si(i 



Si(i 



Ef 1 



E 5 2 (i 
,2 



i=l 

oo 



S 2 (i 



E^ 



E 



S2,i(i 



E 

i=i 

oo 

E 

2 = 1 
OO 

E 

1=1 

E 5 3 (a 
,■2 



E 

j=l 



i=l 

Sx(i)S 2 (i 



S 2 ,i(oo) = 2C 3 , 



53,1(00) = -C4, 



54,1(00) = 3C5 - C2C3, 



5 2 , 2 (oo) = -C4, 



53,2(00) = 3C2C3 ~ 2^5, 



52,2,l(oo) = 2( 5 , 



17 1 



C2C3 + 10C5, 



2 Cs — C2C3, 



11 



<5 - 2C2C3, 



Cs + C2C3 



(72) 
(73) 
(74) 
(75) 
(76) 
(77) 
(78) 
(79) 
(80) 
(81) 
(82) 
(83) 



APPENDIX C: LIST OF MASTER INTEGRALS 



In this Appendix we present the list of master integrals we have used in our calculation. The list updates the result 
that can be found in p0[. The Euclidean master integrals are defined as (for a pictorial representation see Fig. 2): 



h -- 
h -- 

l A - 

h -- 
h = 

ho 
I12 



(-l)p-k 2 \ 
D x D 2 D 3 D 4 D h D^D 7 D 8 / 
1 

D 1 D 2 D 3 D 4 D 7 D 8 
1 

D 2 D 3 D 4 D 6 D 7 D 8 
1 

D 1 D 3 D 5 D 6 D 7 D 8 
1 



1 

D 2 D 4 D 6 D 7 D 8 
1 \ 



D x D 2 D 4 D h r 



h = 
h = 
h = 
h = 
hi 
ha 



1 



D&zDiDnDeDa 
1 

D 1 D 3 D 4 D 5 D 7 D 8 
1 

D 2 D 4 D 5 D 6 D 7 D 8 
1 

D 2 D 3 D 5 D 6 D 7 
1 



D X D 3 D^D 7 

{d 3 d 5 d 6 d 7 ) 
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/l4 ^ 2 5 3 -Da)' 715 \d 3 d 4 d 7 d 8 



ZVWV \£>iD 4 I>5/' 
Il8 = { D 1 D 2 D 5 D 6 D 7 D 8 D 9 ) ' (84) 



where 



(...) = y d D hd D k 2 d D k 3 , (85) 

the momentum p is always considered as incoming and 





= (P- 


kif + l, 


£>2 


= (p - h 


-k 2 ) 2 + l, 


D 3 


= (P- 


ki-k 2 - k 3 ) 2 + 1, 


£>4 


= ip- k 2 


-k 3 ) 2 + l, 


D 5 


= (P- 


*s) 2 + l, 


£>6 


= k l , 




D 7 


- A- 2 






- k 2 




D 9 


= (h- 


-fci-fc 2 ) 2 . 









(86) 



where p 2 = —1. 

/i=C(e) (-5C 5 + ^ 2 C 3 + 0(e) 



,C3 , 1 2 13 4 , J^_2, , , ,_2,__ 4 _ 2 13. 



J 2 = C(e) (^2- + 10C 3 - ^ - -tt 4 + e ^-tt^s + 24C 3 - yCs + 4tt 2 log 2 - -tt 2 - -tt 4 ) + 0(e 2 ) 
T ^ / 1 7 31 103 4, 2 4 /235 28, 2 32, „ o , 8 2 , . , , 

/4= C(e) ( 2 | + 2C 3 + ^ 2 - ^ + e (-^ 12C 3 + 44C 5 + yV - IV) + 0( £ 2 ) 
r 1 3 1 /^ 55 1 2 \ 95 14, 7 2 4 4 

+ e + f Ca. 2 - 42C 3 - f C , + 4. 2 log 2 - §. 2 - §. 4 ) + (e 2 )) , 

T n / 1 7 31 103 2, 1 2 4 4 /235 2 , , 20, „, , 2 3 , , , , 

^ C(e) (^ + 3^ + 3^ + - + 3 C3 + r -45^ +£ br + 3 C ^ 2 + yCa-2< 5 + 4 7r | - 

r ^ n ( 1 3 1/55 1 2 \ 95 8, „ 2 1 4 



+e + 6C37T 2 - 14C 3 - 64C 5 - y> - g^ 4 ) + G(e 2 )) , 



r x / 1 16 16 _ „, 82 /364 200, „„ 2l „ „ , 3 

h= C(e) - ^ - T - 20 + 2 C 3 - -. 2 + e (- - -< 3 + 16. 2 log 2 - 28. 2 - 

+e 2 (1244 + 21C37r 2 - 776C 3 - 126C 5 + 168tt 2 log 2 - ^tt 2 log 2 2 - 188tt 2 + ^tt 4 - ^ log 4 2 - 512a 4 ) + 0(e 3 ) ) 
\ 3 15 3 J J 

T s~i 1 \ ( 2 10 1/26 1 ,\ „ 16, 11 2 /398 248, „ 2l „ 73 2 13 , 

+e 2 (l038 - \c, 3 v 2 - - 96C 5 + 160tt 2 log 2 - ±fn 2 log 2 2 - 129^ 2 + V - | log 4 2 - 512a 4 ) + 0(e 3 ) ) , 

„. .( 1 5 1 / 2 2 \ 10 26, 7 2 /302 94, 2 35 4 

s 2 (-734 + - ^ + 20C3 + f^- 4 + 462C 5 ) + G(e 3 )) , 



21 



Ii2= C(e) (_ + _ + _ + _ +e [-— + — 



280C 3 - f ir 2 log 2 2-^ + 1 log 4 2 + 256a 4 ) + 0(s^ , 



7 253 2501 / 59437 64 2 \ 2 / 2831381 1792,. 256 2 , „ 2272 2 N 
/ll=C(£) b + 2? + 3fe + W + <T29¥-T 7r J +£ \~ 7776 9~ ^ 3~ 2T - , 

, / 117529021 63616 , 9088 2l „ 3584 2l ,„ 49840 2 2752 4 1024, 4 „ 8192 . 

+£ (^6656 27^ 3 + IT* l0g 2 " —** l0g 2 " "81 - + 135-" * " ~ V 2 " — * ) + Oi : ^ 

23 35 275 ( 189 112 
3^2 + Ye + Tsf 

14917 

48 

s A 1 7 25 5 8, / 959 28, 2 

+ ,(.- + fca + 18f5 _^) +0(£3) ), 

/u-c(.,(^ + f + ^ + ^ + ^ + .(-f + »G-^ + ^ 

+e 2 / 14917 + 21QC3 _ 2 + 16 ^ 2 bg2 2 + 145^ 2 _ 62^4 + 8 log 4 2 + ig2 \ + 0(£3 . 

\ 64 3 45 / 

/ 1 7 1 /25 1 2 \ 5 ^ 7 2 / 959 „, 25 2 16 4 

_ 2 / 10493 8 2 5 2 56 4 



~ £ < ~-64~ + 3 C3?r + 25C3 + ?2C5 _ 24^ + 45^ ' + ° (£ } 



/ t . 1 5 1 / 11 1 2 \ 23 16. 5 2 /135 80 A 11 2 37 4 

,,'949 32 2a 176, „„„„ 23 2 37 4 \ ^, 
-e 2 ( — ~ y^Cs - — C 3 - 208C 5 + y^ 2 - -g-TT 4 J + G(e 3 ) 



7 17 = (7(e) ( ™ - ^ - - - 10 - 15e- 21e 2 - 28e 3 + 0{e 4 [ 



e 3 e 2 e 



hs= (7(e) (27r 2 C 3 - 5C 5 + 0{e)) . 

Finally, we should note that integrals I12 and 7i 4 are known to a higher order in e than we have written here; the 
corresponding higher order terms in e can be found in Ref. Furthermore, integrals /13, /15, I\q and 7 17 can be 
performed exactly in terms of Euler gamma functions, as one can see for instance by combining the exact forms in 
Appendix A. This means that these four integrals are known in principle to an arbitrary high order in e. 
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FIG. 2. Pictorial representation of the complete set of primitive 3-loop on-shell integrals I\ - lis. Thick lines are used to 
indicate massive scalar propagators and thin lines indicate massless propagators. 
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